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The effects of the initial temperature in the out of equilibrium quantum field dynamics in the 
presence of an homogeneous external field are investigated. We consider an initial thermal state of 
temperature T for a constant external field J. A subsequent sign flip of the external field, J —> — J, 
gives rise to an out of equilibrium nonperturbative quantum field dynamics. The dynamics is studied 
here for the symmetry broken A($ 2 ) 2 scalar N component field theory in the large N limit. We find 
a dynamical effective potential for the expectation value that helps to understand the dynamics. 

The dynamics presents two regimes defined by the presence or absence of a temporal trapping 
close to the metastable equilibrium position of the potential. The two regimes are separated by a 
critical value of the external field that depends on the initial temperature. The temporal trapping is 
shorter for larger initial temperatures or larger external fields. Parametric resonances and spinodal 
instabilities amplify the quantum fluctuations in the field components transverse to the external 
field. When there is a temporal trapping this is the main mechanism that allows the system to 
escape from the metastable state for large N. Subsequently backreaction stops the growth of the 
quantum fluctuations and the system enters a quasiperiodic regime. 
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I. INTRODUCTION 

Out of equilibrium dense concentrations of particles are present in several important physical systems, as the 
ultrarelativistic heavy ion collisions Q and the early universe These out of equilibrium dense concentrations of 
particles imply the need of out of equilibrium nonperturbative quantum field theory methods, as the large N limit. 

We study here the effects of the initial temperature in the out of equilibrium dynamics induced by a change in the 
external field. We consider the O(N) A$ 4 theory with spontaneously broken symmetry. Initially, an homogeneous 
external field J breaks the vacuum degeneracy and the system is in a thermal equilibrium state at a given temperature 
T. Subsequently, the external field J flips the sign inducing an out of equilibrium dynamics in the system. We study 
this out of equilibrium dynamics using the large N limit method. We pay a particular attention to the effects of 
the initial temperature in the dynamics of the system. This system for the particular case of zero temperature has 
been studied in Ref. [§]. We show here how those results are extended or modified for nonzero initial temperatures. 
On the other hand, the effects of uniform external fields in A<1> 4 theory with broken symmetry have been studied 
in a different framework through classical evolution with random initial conditions Q. Here we study the out of 
equilibrium quantum field dynamics of an initial thermal state after a change in the external field. 

In Section[n|we present the A(<f> 2 ) 2 model and its out of equilibrium evolution equations in the large N limit. The 
evolution equation for the expectation value can be restated in terms of an effective dynamical potential as it is shown 
in Section cm This effective dynamical potential and its maxima and minima help to understand the dynamics. The 
early time dynamics is presented in Section cm where we show that there are two dynamical regimes defined by 
whether the system is temporally trapped close to a metastable state or not. The dynamical regime is determined by 
the values of the initial temperature and the external field, because the two regimes are separated by a critical value 
of the external field, J C (T), that depends on the initial temperature. In this section we also compute the trapping 
time finding that it is shorter for larger initial temperatures or larger external fields. In the next section, Section m 
we study the intermediate time dynamics and find a quasiperiodic evolution with a clear separation of slow and fast 
variables. Finally, we present the conclusions and an appendix where the analytic expression for the spinodal time 
(the trapping time) is derived as a function of the external field and the initial temperature. 


II. THE MODEL 

We consider N scalar fields, $, with a A(<f> 2 ) 2 selfinteraction in the presence of an external field J. The action and 
the lagrangian density are given by 

S = J d 4 xC , (1) 

C = l[d^{x)] 2 - ^m 2 $ 2 --^{$ 2 ) 2 ~ (2) 

We restrict ourselves to the case where the symmetry is spontaneously broken, i.e., m 2 < 0; and we mainly consider 
small coupling constants A, because this slows the dynamics and allows a better study of its different parts. 

We consider here the evolution of a initial thermal state of temperature T after a flip in the homogeneous external 
field J —> — J. A thermal state of temperature T has translational invariance. This translational invariance of the 
initial state is preserved by the evolution equations when the external field J is homogeneous, as in our case. Thus, 
the expectation values of the fields, in particular (<f>) and (<f> 2 ), are independent of the spatial coordinates and only 
depend on time. We want to study the dynamics after the change J —> — J ; this implies that the direction of the 
external field is fixed, and we can choose the axes in the IV-dimensional internal space such that 

f= ((VNJ, 0,...,0) for f < 0, 

\ {-y/NJ, 0,..., 0) for t > 0. 

For an initial thermal state we have an expectation value parallel to the external field. The O(N) invariance of the 
lagrangian © for J = 0 together with the fixed direction of J guarantees that the expectation value remains parallel 
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to the external field during the evolution. Therefore, the following decomposition can be done 

$(ir) = (cr(x),Tv(x)) = (VN<j>{t) +x(x), i ?(x)) , (4) 

with \/N(/)(t) = (a(x))', thus, (x(x)) = 0. While in the remaining TV — 1 directions transversal to the expectation 
value, (7?(x)) = 0. 


A. Evolution equations in the large TV limit 

We present here the main concepts and results that leads to the derivation of the evolution equations in the large 
TV limit. More details can be found in Refs. Hi 

As we have one direction parallel to the expectation value and TV — 1 transversal, the fluctuations in the transverse 
directions dominate in the large TV limit, while those in the longitudinal direction only contribute to the evolution 
equations as corrections of order 1/TV. 

In the Heisenberg picture we have 


7r (x, t) 


d 3 k 

V / 2(2tt) 3 




+ 


— ik-x 


(5) 


where a^, at are the annihilation and creation operators respectively, that satisfy the usual canonical conmutation 
rules. The functions (p%(t) are the mode functions of the ff field. 

The evolution equations for t > 0 are 


4>(t) + Ml(t)<j>(t) = -J , 

0jg(*) +vl(t)<p u {t) = 0 , 


where we define the effective frequency wg as 


w 2 u (t) = k 2 +M 2 d (t), 


( 6 ) 

(7) 

( 8 ) 


Mi = to 2 + 


4>\t) 


(* 2 m 


TV 


(9) 


The last term in the effective mass squared, A! 2 , gives the quantum and thermal contributions. It can be interpreted 
as the mean effect due to the 7 ? particles, and it is given by 


<^> 

TV 


d 3 k 


S d = -- 


J M) 3 

9{k — k) 
2 k 3 


M 2 d {t) , 


z coth (^M 


-Sd 


( 10 ) 

( 11 ) 


with M 1 = ksT , and k an arbitrary renormalization scale that we choose, for simplicity, equal to the renormalized 
mass, k = \mn\. (For details on the renormalization procedure, that leads to the subtraction S d see Ref. §5].) 
The previous equations are already written for the renormalized magnitudes. It is important to stress that the 
renormalization is independent of the temperature {?]. The hyperbolic cotangent factor in the Eq. m is due to the 
initial thermal state (see for example Ref. 0]). In addition, the fact that the initial state is a thermal state implies 
the following initial conditions for the expectation values and the modes: 


m = 00 , 0(0) = 0 , (12) 

V’g(O) = J-j-r » 0fe(O) = ~Wuk{ 0), (13) 

VWfc(0) 

where 0o is the stationary solution for t < 0 with minimal total energy. The evolution equation for t < 0 is the Eq. © 
changing J for — J, therefore the stationarity condition for t < 0 is M 2 d (j >o = J. Note that the initial conditions for 
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the modes are spherically symmetric in the momentum space and the evolution keeps this symmetry (because the 
evolution equations are also spherically symmetric). 

We introduce the following adimensional variables: 


r = \m\t, 


J 


A J 


k A Md(t) 

q = I—r i 9 = n ] • A/f ( r = | | 

m o7r z m 


2N \m\ 3 ’ 


- 5 S ( r ) = Ap AA ’ 


5 = 




in terms of these adimensional variables the evolution equations for r > 0 are 

v(t) + = -j , (p q [r) + w 2 (r)(p q (r) = 0 , 

where 

M 2 (t) = -1 + t? 2 (t) + pS(r), w 2 = q 2 + M 2 (t) , 


3£(t) = g j q‘ 


dq 


\ L Pq{ T )\ 2 coth 




-5 


q 2q A 

The initial conditions O and are expressed as 

V(0) = do , g(0) = 0 , 


<Pg(0) = 




<Pq( 0) = -iyjuJq{Q) 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


III. EFFECTIVE DYNAMICAL POTENTIAL FOR THE EXPECTATION VALUE 


We can define an (adimensionalized) energy for positive times as the (adimensionalized) expectation value of the 
T 00 component of the energy-momentum tensor. 


A 


e r = 


ry 4 ry 2 , 1 2 


2W |„F (T ) = T-f + 2^-T 


gZ , ( 9 S)= 


-jq- 


■ 9 Hl 2 »ud^« (0) 


g 


+ k 


| gjq | 2 coth 


q 

PVq{ 0) 


-Si 


(23) 


-5 2 


where 5 l and 52 are two subtractions due to the renormalization 


51 = q + 

5 2 = - — 


A4 2 6(q — 1) 


2g 8g 3 
1 AW 2 6>(q - 1) 

q 2q 3 8 q 5 


(m 2 ) 2 + 

3 (M 2 ) 2 


d 2 M 2 
dt 2 

d 2 M 2 


dt 2 


(24) 


For a time independent external field j, as it is the case for r > 0, the energy is conserved. From the expression iPl) 
we have in the (rj,gY,) plane the restriction e > V de -,T>o{jh £), where 

^ i 

(25) 


T , / _ V 4 V 2 , 1 2 v 5 s , (ff s ) 2 , 1 , ■ 

V de -r>o{v, Y>) = — - — + -g gT,- — + —— + - + J?y ■ 


^ / de;r>o( ? ?) £) can be interpreted as a dynamical effective potential because the evolution equation 11711 for ry can be 
written as 

d 

V( T ) = --Q^V de . T> o (ry,£) . 

It must be stressed that V de;T> o is an effective potential only for g (and not for the modes). 


( 26 ) 
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A. Equilibrium states for the effective dynamical potential 

The dynamical effective potential for r < 0 can be defined as Vd e - T <o{jh^,) = Vde;r>o(v, E) — 2 jg- Therefore, the 
stationary states for the initial dynamics for times r < 0 (before the external field has been flipped) are the solutions 
°f Vde-,T< o (rj) = 0 (the apostrophe means g derivative). Thus, they verify 


rf + (~l + gT,)rj- j = 0 . 


(27) 


For small external fields. 


3 < 3d = 2 


( 1 -gS) 
27 ’ 


(28) 


we have three roots. There is a global minimum, that corresponds to a stable equilibrium state, at the value of r) 

j 


n 1-gT, 1 

Vo = 2y ——— cos < - arccos 


27 


2y (1 -9X) 3 

= (1 - gS) 1 / 2 + 1(1 - gT.)~ l j - 1(1 - g^)- 5/2 j 2 
+ 1(1 - gT,)~ 4 j 3 - 111(1 - ffS)- 11 / 2 / + O(f) , 


(29) 


a local minimum (metastable equilibrium state) at 


'1 -gT, 1 

Vi = 2 \/ —-— cos < - arccos 


27 


(1 - gS)3 


27r 

T 


= -(1 - &) 1/2 + 1(1 - gS)- x j + §(1 - gX)~ 5/2 j 2 


(30) 


+ 1(1 - 5 E )- 4 j 3 + 1|(1 - 5 E )- 11 / 2 i 4 + 0(j5) J 


and a local maximum (unstable equilibrium) at 


„ /1 — oE 
V 2 = 2 W —-— cos ■ 


■arccos 


27 


2 V (1-gS) 3 


47T 

T 


(31) 


= -(1 - g£) _1 j - (1 - gS)- 4 j 3 + o(j 5 ). 


On the other hand, for large external fields (j > j^) there is a single extreme, that corresponds to a global minimum 
(stable equilibrium), at the value 


go 


-'0 


1/3 


1 + f/l- 


4(1 ~~ ffE ) 3 
27f 


1/3 


+ 

N 


4(1 - 5 E ) 3 


27j 5 


1/3 


(32) 


= A 3 + l(i - »s)r 1/3 - lj(i - s£) 3 r 5 / 3 +o<j- 7/3 ). 


We consider here the more interesting case, j < jd , where a potential barrier is present. The initial conditions are 
those of the stable equilibrium state of the dynamics for times r < 0. In Fig. ^ we show the location of the minima 
and the maximum of Vd e -,r<o(v) for gE(0) = 0.05 and an external field j = 0.2 < jd- 
For zero initial temperature </E(0) ~ g <C 1. Thus, the zero temperature results of Ref. 0 are recovered neglecting 
gYi in Eqs. E3 03 

IV. EARLY TIME DYNAMICS AND DYNAMICAL REGIMES 

The fact that the initial state is not the ground state but a thermal state increases the initial value of g E. It can 
be shown that E(0) is approximately given by 


s(o) = e t= 0 (o) + — /r 2 

o 


( 33 ) 
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FIG. 1: Dynamical effective potential V*, ;r < o as a function of 77 for g£(0) = 0.05. 

The value of the external field is j = 0.2. The positions of the global minimum 770 , 

the local minimum rji (metastable state), and the relative maximum 772 are shown. 



FIG. 2 : Dynamical effective potential Vd e -, T < 0 and Vd e -, T >o as a function of 77 for <?E(0) = 0.05. The value of the 
external field is j = 0.2. The positions of the global minimum r/o , the local minimum 771 , and the local maxi¬ 
mum 1 J 2 are shown. There is a potential barrier (j < jd), and the systems gets temporally trapped ( j < j c ). 


£ T_o (0) is the zero temperature value, that for j <C 1 only depends on j (see Table 0. The second term on the 
righthand side of Eq. m is the thermal contribution computed in the hard thermal loop approximation (O'], i.e., 
assuming the main contribution comes from the modes with q ~ /3 _1 and we are in the case /3 _1 Ad 2 (0). [On the 

other hand, for /? _1 -C Ai 2 (0) thermal effects are negligeable, because /3 -1 <C Lu q and coth(/3w g /2) = 1 + 0(e~ !)ulq ).] 


j 

0.01 

0.05 

0.10 

0.20 

0.25 

- 10 2 E' J=o (0) 

1.2 

4.2 

6.5 

9.6 

11 


TABLE I: £ T-O (0) values obtained for various values of the external field j (j C 1). 

After the initial flip of the external field sign at r = 0 the positions of the absolute minimum and the relative 
minimum of the potential are interchanged, and the state of the system becomes a metastable state (Figs. 0 and 0. 
The height of the potential barrier is Vd e -,T>oW 2 )-> where rj' 2 = —772 is the maximum of Vd e -,T>o(v)- The existence of 
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FIG. 3: Dynamical effective potentials Vd e -,r<o and Vd e - T >o as a function of 77 for <?E(0) = 0 . 2 . The value of the external 
field is j = 0.2. The positions of the global minimum 770 , the local minimum 771 and the local maximum 772 are shown. 
The potential barrier is not high enough (j < j c ), and the system rapidly evolves towards the absolute minimum. 


this potential barrier gives rise to two different dynamical regimes. In the first one the system can directly overcome 
the barrier Vd e - T >o(vo) > Vde-r> 0 (^ 2 ), and rapidly reaches the neighborhoods of the global minimum. While in the 
second regime the system can not overcome the barrier directly Vd e -,r>o(vo) < Kze;r>o(?? 2 ): and it gets temporally 
trapped close to the metastable state. Solving the equation 

Vde;T>o{Vo) = Vde\T> o(%) (34) 

the external field critical value j c that separates the two regimes, untrapped \j\ > j c , or trapped \j\ < j c , is obtained 


■t= o 


Jc = 


(1- 5 E(0 )) 3 / 2 


with 


Jc T=0 = \ 2 


(l3 2 + 15^5) 


19 3 


= 0.243019... ; 


(35) 


where we have used that for very early times we can make the approximation gTi(r) ss gE(0). An explicit analytic 
expression for the critical external field as a function of the initial temperature is obtained using Eq. iTOll . 


3c(P) 


•T—0 



with 



(36) 


This analytical formula is compared with the numerical results for g = 10 -2 and g = 10 ~ 3 in Fig.QJfinding a good fit. 
The two dynamical regimes (untrapped \j\ > j c and trapped |j| < j c ) correspond to the two regions in the (j, (3~ l ) 
plane separated by the j c (P) curve. This curve can be better understood recalling the effects of the temperature 
and of the external field on the potential. As the temperature increases gT, grows, and the barrier in the 77 direction 
between the two minima of the potential diminishes until it disappears. The barrier disappears even for j = 0 when 
fj~ l > (/3^ =0 ) _1 . On the other hand, the effect of the external field j is to tilt the potential; for j > jj =0 the potential 
is so tilted that there is no temporal trapping even for (3~ l = 0. Therefore, Eq. states how the trapping can 
disappear due to a combination of both effects. 

We analyze now the two dynamical regimes: 


A. j < j c {/3) 

In this regime (Fig. [21) the system does not have enough energy to jump over the potential barrier and it remains 
some time oscillating around the metastable minimum (Fig. 0. The system is only temporally trapped because the 
spinodal instability of the modes (due to the fact that A4 2 (r) < 0) makes the amplitude of the modes grow. This 
growth allows the field to surround the maximum of the potential in the N — 1 directions transversal to the external 















FIG. 4: Critical external field j c as a function of /3 -1 for g = 10 -2 and g = 10 ~ 3 obtained from numeri¬ 
cal simulations (dots) and from the analytical formula in Eq. 1361 1 (full line). There are two regions in the 
(j,/3 _1 ) plane corresponding to the two dynamical regimes: for j < j c the system is temporally trapped in the 
metastable state; for j > j c the system rapidly reaches the neighborhood of the global minimum of the potential. 



T 


FIG. 5: i) as a function of time for j < j c (/3). It corresponds to the values g = 10 3 , j = 0.2, (3 1 = 2. 


field, and finally the system reaches the neighborhood of the stable minimum. (Another effect that leads the system 
to the global minimum is tunnelling. However, our computation neglects this effect because it only takes place in one 
of the N internal directions and thereby it is subleading in the large N limit.) 

The classical evolution equations (recovered for E = 0) predict that the expectation value rj oscillates forever 
between the value 770 and the turning point value rj r given by 


brfe;r>0 ijlr ; ^ — 0 ) — Fde;r>0(^?05 ^ — 0 )- 


(37) 


This is also the case when gE(0) <C 1 (otherwise the trapping is very short). Solving Eq. (13711 order by order gives 




11 


7 -f 


1049 

128 


j 4 + 0(j 5 ). 


(38) 


In Fig. M we plot the behavior of the squared effective mass. For early times, M 2 oscillates with a negative 
average value. This makes the low momentum modes spinodally unstable and they grow exponentially. This implies 
a quasiexponential growth of g E for early times, and finally gT, becomes of order one (see Fig. □) This makes that 77 
is no longer trapped close to the metastable minimum (see Fig. 0 . 
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FIG. 6: AI 2 as a function of time for j < j c (/3). It corresponds to the values g = 10 3 , j = 0.2, (3 1 =2. 





FIG. 7: gT as a function of time for j < j c {(3). It corresponds to the values g = 10 3 , j = 0.2, (3 1 = 2. 


Parametric resonances are also present on the evolution, due to the oscillations of AI 2 , implying the transfer 
of energy to certain momenta that grow in amplitude. This two mechanisms, spinodal instability and parametric 
resonance, transfer energy from the expectation value rj to the modes tp q increasing gT,. However in this case, \j\ < j c , 
the main contribution is the spinodal instability. 

A detailed analysis of the spinodal instability allows to obtain the trapping time, or spinodal time r s . (Here we 
summarize the results, the details can be found in Appendix fXIl At early times (r < r a ), the effective squared mass, 
AI 2 , oscillate with a negative average. Thus, we can define an effective average squared mass —/j, 2 , and in the modes 
evolution equation replace AI 2 (r) by —/i 2 to obtain an estimate for \ip q \. This is a good approximation after 
several complete oscillations of AI 2 (r) (i.e., r > 1/g) 1 and when the time dependence of AI 2 (r) is slow (adiabatic 
approximation). The effective average squared mass is given by [Eq. (IA4I) 1 

= + ( 39 ) 

The quasiexponential growth of gT(r) for r < t s is 

J JJL9^ = °( T ) for /3 _1 » V > 

l ffS^=°(r) for (3- 1 < /j, , 


9^s(r) 


(40) 
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a 

9 

2-10 2 

2-10 3 

2-10 4 

2-10 5 

2-10 6 

2-10 7 

t s (numerical) 

4.8 

9.6 

11.5 

15.1 

16.7 

20.1 

Ts (analytical) 

7.1 

10.1 

12.8 

15.6 

18.2 

20.8 

error (%) 

33.1 

4.1 

10.8 

3.0 

8.4 

3.6 


TABLE II: Spinodal time values for different values of — with fixed f3 = 0.2 and j = 0.2. The values has been obtained 
from the numerical simulation of the complete evolution equations, and from the analytical approximation Eq. 113 - 


wit h 




g^Jle 2 ^ 

8r 3 / 2 


(41) 


the value for zero temperature. 

After a certain time, the spinodal time r s , the quantum and thermal effects start to be important in the dynamics, 
gE s (r s ) compensates —g? and the exponential growth of the mode functions stops, then the mode functions start to 
have an oscillatory behavior. Thus, the spinodal time r s is defined as 

gY> s {r s ) = g 2 . (42) 


A good approximate expression for the spinodal time is 


T s 


b [ lo s (ife) - lo e fe) +1 lo s(^) 


for [3 1 > n , 
for /3" 1 <C n . 


(43) 


We see that the spinodal time for /3 _1 <C g is the same as that for zero temperature rj =0 (see Ref. jj|). On the other 
hand, for /? _1 /i we see that increasing the initial temperature the spinodal time decreases. Therefore, a higher 
initial temperature implies a shorter trapping period, that is consistent with the fact that for greater temperatures the 
initial transversal fluctuations are larger and the system is closer to the end of the trapping period. The contribution 
of the transversal fields grows rapidly following Eqs. Efli-EB. turning around the maximum of the potential and 
reaching the neighborhood of the absolute minimum at earlier times. Increasing the external field j also shortens the 
trapping period, as can be shown in Eq. El (M ~ Vj)- 

In Table EH we compare the numerical results obtained solving the evolution equations © U3| and applying the 
t s definition in Eq. (EH . with the values predicted by the analytical expression in Eq. El- We see a good agreement 
(with discrepancies of 10% or smaller) except for low ^ values (that correspond to high temperatures). This later 
discrepancy is not troublesome because the expression El has been obtained for times r > F and this condition does 
not hold for high temperatures, as the the spinodal time decreases for decreasing In fact, when j = 0.2 we have 
g 2 ft 0.23 [Eq. (1-391) 1 and therefore F ft: 2.1; as for ^ = 2 • 10 2 we have a spinodal time r| lm ft 4.8 ~ 2.1 fTablelTTIl. 
Eq. 11431) can no longer be applied. However, for greater spinodal times (that corresponds to ^ > 2 • 10 3 ) the behavior 
predicted by the Eq. EH is reproduced with a good agreement, as we show in Table ITT! 


B- j > jc(P) 

The main characteristic of this regime (Fig. © is that the system has enough energy to overcome the potential 
barrier directly. Therefore no trapping is present in this case. 

Also in this case there is energy transfer between the expectation value and the modes of the field. The low 
momenta modes are spinodally unstable and grow exponentially only during the time intervals when A4 2 (r) < 0 
(Fig.©. Parametric resonances are also present due to the oscillations of the effective squared mass. The Figs. ©© 
and d reflect the system dynamics in this regime. Both mechanisms transfer energy from the expectation value rj to 
the quantum fluctuation ifk- Therefore, ffE(r) increases while the amplitude of the oscillations of g decreases. 













0 5 10 15 20 25 30 35 40 45 50 

x 

FIG. 8 : 77 as a function of time for j > j c (/3). It corresponds to the values g = 10 -3 , j = 0.2, /3~ 



0 5 10 15 20 25 30 35 40 45 50 

X 

FIG. 9: 77 as a function of time for j > j c (/3). It coorresponds to the values g = 10 -3 , j = 0.2, (3~ 



0 5 10 15 20 25 30 35 40 45 50 

x 

FIG. 10: gT, as a function of time for j > j c ((3). It corresponds to the values g = 10 -3 , j = 0.2, /3 
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V. INTERMEDIATE TIME DYNAMICS 

After the early period described in the previous section, the system enters in a quasiperiodic regime. This behavior 
has been already observed for the particular case T = 0 in Ref. y], while for different initial conditions the oscillations 
damped much faster HOI 

This quasiperiodic behaviour found for t](t) and £/E(t) suggests that these quantities are approximately governed 
by an effective Hamiltonian with a few degrees of freedom. Actually, we found (as in Ref. 0 ) that the full numerical 
solution of Eqs. m-m gives that t](t) and sE(r) are approximately related as 

gE(ri) = 1 + c 0 j - ( l-c 1 j){rj + jf (44) 

where Co and C\ are positive numbers of the order j° and g° for small g, that are obtained fitting the numerical 
solution. See Fig. and Tables ITTT1 and CD 


j 

Vgr 1 

0.05 0.10 0.20 

0 

0.56 0.54 0.40 

0.1 

0.48 0.48 0.37 

0.2 

0.42 0.43 0.20 

0.3 

0.40 0.40 0.11 


TABLE 

III: 

co 

values 

obtained fitting gE( 77 ) 

to Eq. 


in the 

time 

interval 

r G [800, 1000 ] 

for g 


10~ 2 

for 

various values of the 

external 

field 

j and 

the 

initial 

temperature. 





j 

0.05 0.10 

0.20 









s/gp- 1 











0 

1.49 1.02 

0.55 









0.1 

1.62 1.12 

0.63 









0.2 

2.03 1.42 

0.98 









0.3 

2.91 2.04 

1.49 





TABLE 

IV: 

Cl 

values 

obtained fitting gE,(g) 

to Eq. 

iTHli 

in the 

time 

interval 

r G [800, 1000] 

for g 

~ 

10” 2 

for 

various values of the 

external 

field 

j and 

the 

initial 

temperature. 


Using this approximate expression the effective squared mass can be expressed as 

A4 2 (ry) = -1 + if + gE(g) = j [c 0 - j + aj 2 - 2(1 - cij) g + ci rf] . 
The evolution equations ® nu in this approximation take the form, 

ij + M 2 (g)g = -j . 

We find integrating on 77 , 

1 


where, 


:V T j Vint (v) — j Eint 


V int (v) = 11 + 77 (co — j + c\j 2 )rf - ^(1 - Cl j) if + J? 7 4 , 


Ejnt — Vj n i ( 771 ) • 


(45) 

(46) 

(47) 

(48) 

(49) 


Notice that Ei nt depends on the initial conditions, and that 771 is a turning point of the motion. Eq. (IT71) can be 
integrated as follows, 


\/27( T -n) = 


dr] 


?7i y/^int Vint (^) 


(50) 
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FIG. 11: Trajectory in the ( 77 , 5 s) plane for r < 1000 and 

[800, 1000]) fit to Eq. (151 . corresponding to the values g = 10 -2 , 


intermediate time (r 
j = 0.20, (3~ l = 


e 

1 . 


with rji = 77 (ri). 

The fourth order polynomial Ei nt — Vi n t{rf) has always two real roots 771 < 772 corresponding to the turning points. 
Depending on the value of j, the two other roots can be: 

• i) a pair of complex conjugate roots t/r ± 777 /, for which we define 


VR - m , = Vi 

(VR - Vi ) 2 + V 2 ’ ~ (VR ~ Vi ) 2 + V 2 ’ 


(51) 


• ii) or a pair of real roots 771 < 772 < 773 < 774 , we then define: 


2 _,2 = 4(773 - 772)(774 - 773) 

V3 ~ Vi Vi-rn' ~ (V 2 ~ Vi){V4 ~ Vi)(V3 ~ Vi) 2 ’ 


(52) 


We also introduce two other quantities to simplify the formulas: 

d=X = [(a — d) 2 + b 2 ] 1/4 ; C = X ^2\M 2 ( Vl ) Vi + j\ (53) 

772-771 

It is convenient to use u = v l m as integration variable in Eq. lt50l) . The solution of Eq. lf50l) can be expressed after 
calculation as 


, , ,_ (772 - 77!) [1 - cn(C l (T - Tj), k)} _ 

77 T 11 1 + (772 - ? 7 i)* 2 - [1 - (772 - 771 )X 2 ] cn (C(t - n), k) 


(54) 


where cn(^, k) is the Jacobi cosine function, and 


k 



a — d 


the elliptic modulus. 

The solution (T >1 oscillates between 771 and 772 with period 


T=l/T(fc) 


(55) 


(56) 


where K{k) stands for the complete elliptic integral of the first kind and C is given by Eq. (15311 . The analytical and 
the numerical solutions are compared in Fig. 1121 
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FIG. 12: Quantum evolution for the expectation value g{r) (full line) compared to the analytical intermediate time 
approximation of Eq. (ICT (dashed line). Corresponding to g = 10 —2 , j = 0.20 and /3 _1 = 1 [for this case j < j c (P)}- 


T 

400 700 1000 

j Ei n t {j = 0.05, vs/3 _i = 0.0) 
jEintij = 0.05, Vs/3 -1 = 0.1) 
jEintU = 0.05, Vs/3 -1 = 0.2) 
jEintij = 0.05, Vs/3 -1 = 0.3) 

0.022 0.016 0.015 

0.019 0.013 0.009 

0.015 0.011 0.006 

0.012 0.007 0.006 

jEintij = 0.20, Vs/r 1 = 0.0) 
jEintij = 0.20, Vs/?- 1 = 0.1) 
jEintij = 0-20, Vs/3 -1 = 0.2) 
jEintij = 0.20, Vs/!" 1 = 0.3) 

0.011 0.004 -0.011 

-0.004 -0.014 -0.024 

-0.027 -0.035 -0.046 

-0.050 -0.058 -0.064 


TABLE V: 


E;, 


[Eq. dm ] for g = 10 2 and various values of the external field j and the initial temperature. 


As in the zero temperature case (Ref. Q), the numerical solution of the large N limit evolution equations 117M20II 
would be periodic if Ei nt were exactly conserved; however it is slowly decreasing implying a slow damping of the 
oscillations. This damping is smaller the smaller is j. See Table El 

Therefore, we have seen that the system at intermediate times presents a clear separation between fast variables 
and slow variables. g{r), gT(r), and A4 2 (t) oscillate fast, while E int slowly decreases. We have found the fast time 
dependence of the order parameter and the quantum fluctuations in terms of rational functions of Jacobi cosines. 

The analytic expression for gT(r), and A4 2 (t) are obtained from the relations [Eq. (1441) 1 and A4 2 (rj) [Eq. 

(TTsli l. just using the analytical solution for g{r). The relation gTig) indicates that gT and g oscillate with opposite 
phase, up to order j terms; gT{g) = 1 — g 2 + O(j). This O(j) terms are very important because they imply 
A4 2 (t) = — 1 + g 2 + gT, = 0(j). Therefore the time averaged squared mass is positive, and tends asymptotically to a 
positive value. This can be seen also from Eq. ESI . Thus, in the limit of zero external field we recover a zero effective 
mass squared consistently with the presence of Goldstone bosons for zero external field. [Recall that for j = 0 in the 
broken symmetry case the squared mass goes to zero for initial energies below the potential at the local maximum 

(77 = 0 ) 601 .]“ 


VI. CONCLUSIONS 

We have studied here the effects of the initial temperature in the out of equilibrium dynamics induced by a flip in 
the sign of the external field. The study has been done in the leading order of large N for the A($ 2 ) 2 theory. 

Before the flip of the external field, the system is in a thermal state, with the expectation value of the field in the 
minimum of its effective potential. After the flip the system evolves in either of two dynamical regimes. The difference 
between the two regimes is the presence or absence of a temporal trapping of the system in the neighborhood of a 
metastable state. The presence of this trapping close to a metastable state is reflected in the fact that the expectation 
value oscillates without changing sign (despite of the sign change of the external field). 
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From the evolution equation in the large N limit, we have found a dynamical effective potential for the expectation 
value in the direction of the external field that helps to understand the dynamics. The trapped regime corresponds to 
the case where the system cannot overcome the barrier appearing in the dynamical effective potential. We have shown 
that the dynamical regime is determined by the initial temperature and the external field. There is a temperature 
dependent critical value of the external field that separates both regimes, see Eq. EH). For external fields smaller 
than the critical one, the system is temporally trapped close to the metastable equilibrium position of the effective 
potential. For external fields greater than the critical one, the neighborhood of the stable equilibrium is rapidly 
reached and no trapping is observed. 

Spinodal growth of the quantum fluctuations is the main mechanism that makes the trapping be temporal. The 
increase of the quantum fluctuations diminishes the potential barrier between the two minima of the dynamical 
effective potential for the expectation value. We have found an approximate analytic expressions for the trapping 
time (spinodal time), see Eq. IITill . It shows that the trapping time becomes shorter for larger initial temperatures or 
larger external fields. 

After the trapping ends, the transfer of energy from the expectation value to the quantum fluctuations takes place 
through two mechanism: spinodal instabilities and parametric resonances. The spinodal instabilities appear during 
the time intervals when the effective mass is negative, while the parametric resonances are due to the oscillations of 
the expectation value. 

We have found that these processes lead for intermediate times to a quasiperiodic regime with a clear separation of 
slow and fast variables. We have explicitly solved the time evolution for the fast variables that oscillate periodically, 
while the other variables change slowly due to a tiny damping of the oscillations. The effective squared mass oscillates 
around a positive value of the order of magnitude of the external field. Thus, in the limit of zero external field we 
recover a zero effective mass squared consistently with the presence of Goldstone bosons for zero external field jf|. 

An interesting question that remains still open is which of these results would be modified by the next to leading 
order terms of the large N approximation |ll| and how much. The inclusion of next to leading order terms allows 
to include the effects of longitudinal fluctuations. (The longitudinal fluctuations have quantum and thermal origin, 
and later they reflect the effects of thermal jumping over the barrier and of quantum tunnelling through the potential 
barrier.) This inclusion will allow to compute the probability of finding the system close to the stable minimum before 
the spinodal time, and it will also give corrections to the spinodal time; these corrections are expected to be small, as 
the spinodal time will be mainly determined by the quasiexponential growth of the transversal fluctuations. Also the 
inclusion of next to leading order terms have been shown to damp faster the oscillations on the magnitudes during 
evolution EH- Therefore, how much the next to leading order terms will increase the damping rate of the intermediate 
time quasiperiodic behavior is still an open question. 
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APPENDIX A: SPINODAL TIME 


We present here the early time solution for the modes in the spinodally resonant band, and an estimate for the 
spinodal time, r s , in the case where the external field is smaller than the critical value, j < j c (/3). In this case the 
spinodal time corresponds to the time the system is trapped close to the metastable state. 

At early times (before the spinodal time, r s ), the effective squared mass, A4 2 (r), oscillates with a negative average. 
Thus, we can define an effective average squared mass — /r 2 and replace A4 2 (r) by — /i 2 in the modes evolution 
equation El to obtain an estimate for \<p q \. This is a good approximation when the time dependence of AI 2 (t) is 
slow (adiabatic approximation), and after several complete oscillations of A4 2 (r). Assuming <?E(t) <C 1 (otherwise r s 
is very short), we have M. 2 (t) rts — 1 + 77 2 (r). As g evolves according to Eq. (TTTIl with a frequency M(t ), both rj and 
AI 2 (r) oscillate with a period of the same order V Therefore, as the effective average squared mass approximation 
is valid when we average over several periods, this requires 





(Al) 


(We will show at the end of this appendix that requiring t s > 1/g implies gT,(r) -C 1 during most of the time interval 
r G [0, T s ].) 
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The average squared mass, — /i 2 , is estimated as the average value of M 2 between its maximum A4 2 nax , and 
minimum M^ in value. These values correspond to the g field oscillations between g 0 [Eq. (1231) ]. and the turning 
point value rj r [Eq. 13811 ] . Lets call rp the time over which the average is done, it results 

-V = — [ Ad 2 (r) dr ~ i (M 2 min + M 2 max ) ; (A2) 

'P J o z 

Neglecting gY for these times r < rp, 

•WL » -1 + rfi » J ~ \f + \f ~ f + 0(j 5 ), 

1 23 ( A3 ) 

M 2 mm « -1 + g 2 « — 3j - 2/ - yi 3 - 4j 4 + 0(j 5 ). 

Using these expansions on j we get from the Eq. 

V = -3 - \f - \f -\f + 0(f). (A4) 

The evolution equations for the modes are then 

(^2 +f ~ M 2 ) Vq( T ) = 0 ( A5 ) 

with the initial conditions 


<Pq( 0) = (g 2 + Ad 2 (0)) 1/4 , ip q (0) =-i(q 2 + M 2 {0)) 1/4 . 


(A6) 


For low momenta (those in the band 0 < q < g) the effective squared frequency q 2 — g 2 is negative (spinodal instability) 
and the <g q {r ) have an exponential behavior. 

Solving this second order ordinary differential equation, 


<Pqif = 


1 


g Ty/lJ, 2 -q 2 
g— r V 7 M 2 —<? 2 


(A7) 


2 + Ad 2 (0)) 1 / 4 

(vV-9 2 -*vV+AI 2 

+ (\/M 2 - q 2 + if Q 2 + M 2 

The contribution of the spinodal band to gT, is 

gY, s {r)=gj dqq 2 \tp q {r)\ 2 coth ^ j . (A8) 

This spinodal band gives the dominant contribution to gE s (r). The exponential growth of these low momenta modes 
makes that after a short time they give the main contribution. 

Inserting Eq. m in Eq. IA8I) we have 


gY, s (T)=g dq < q 


g 2 + M( 0) 2 


V (l - $) x /g 2 +M(0) 2 


2 TfJ.\ 1- 


coth 


0y/q 2 + M(O) 2 


(A9) 


where we have neglected the exponentially decreasing term of ip q . This is justified by the condition ED- 
It is convenient to separate three dynamical regimes: 

• High temperatures, /3 _1 g/g, this implies gll( 0) > 1, i.e., j > j c (P), and there is no trapping. 

• Medium temperatures, gjg /3 _1 g, this implies gE s (r) _0 (r) as we show later. 

• Low temperatures, /I -1 <C g , this implies gE s (r) ~ gTjJ =0 (r) as we show later. 

Therefore, here we are only interested in the medium and low temperature cases where the system can be trapped 
close to the metastable state giving a nonzero spinodal time. We devote the following subsections to these cases. 
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1. Medium temperatures (g/g /3 1 /§> g) 

We call medium temperatures those that verify /3 _1 uj q (0) for all the modes in the spinodal band, i.e., 


and also give gT,( 0) <C 1, that implies 


/3 1 > /i ; 

/3 _1 <C g/g . 


This allows the following simplifications 

M 2 (0) = -1 + 77 2 (0) + ff E(0) « -1 + 77 2 (0) , 

and also in this regime (/3 _1 g) 


coth 


Pvq( 0) 


= coth 


' py/q 2 + M 2 (0)' 


Py/q 2 + M 2 ( 0) 


Therefore, 


gT, s (T)=g dq\q 


g 2 +M( 0) 2 


V (l - $) y/q 2 +M(0) 2 


2 rgj 1 - 


(A10) 

(All) 

(A12) 

(A13) 

(A14) 


/V* +A4 2 (0) 

In addition, A4 2 (0) is of the same order of magnitude as g 2 , and we make the approximation A4 2 (0) ~ g 2 . We get 

1 


gT, s {r) = ^J^ dqq 


2 Tfl\ 1 - 


M 2 1- 


1 + K 

fl 2 


(A15) 


It can be shown that the greatest contributions to the integral come from the integrand values that have q = O(0.1g) 
(the exponential enhance low momenta while other factors suppress very low momenta). This justifies making a 
Taylor expansion to second order in —, both in the exponential and in the previous factor. The result obtained is 


9^s{t) = Tj / dq —r e V ) 

P Jo V 


(A16) 


and making the change of variable f = q. — the previous expression becomes 

V A 4 


v f N 9 e 
9^s{t) = - 


2 th r\/W 


(3 7 gr 3/2 J o 




2 p -«; 


(A17) 


For gr —> oo the integral gives y / 7r/4, that is a good approximation for /jr > 1 [condition llAHl ], Therefore, a good 
approximate expression for the growth of the quantum and thermal contributions in the spinodal band is 


9^(t) 


9 /7T e 


3 2r/i 


(A18) 


4/3 Y p, r 3 / 2 ’ 

that is valid for times - < r < r s and temperatures g «C /3 _1 <C -. 

A 4 ^9 

In the case of zero temperature we have (/E^ =0 (r) ~ gv ^/2 " (see Ref. 0 ) and we can express the result obtained 


as 


. 


(A19) 


After a certain time, the spinodal time r s , the quantum and thermal effects start to be important in the dynamics, 
//E s (r s ) compensates — g 2 and the exponential growth of the mode functions stops, and they start to all have an 
oscillatory behavior. Thus, the spinodal time r s is defined as 


gT, s (T s ) = g 2 . 


(A20) 

















Using Eqs. in and (IA20I) we have 
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and taking logarithms we obtain 


implying 


g /?e 2T ^ 2 

4 ^ ’ 


T s 




'log {/JT S ) 


T s = — 


( 4 /fy \ 
. 0 (****? 


2/i 


4/x 


log 


■log 


4^\ 

s-A/ 




!°g^ 


(A21) 


(A22) 


(A23) 


this expression is related with the spinodal time for zero temperature, tJ (see Ref. i) through the equation 

t s = rj =0 - log + ° ^og lo g ~ ) ■ (A24) 

One direct conclusion from the expression 1 A 241 is that an increase of the initial temperature diminishes the spinodal 
time. 


2. Low temperatures (/3~ _1 


For temperatures that satisfy 


/3 1 <C g 


we can make the following development for the hyperbolic cotangent, 

'/3w ? (0)' 


coth 


= coth 


0y/q + M 2 (O) 


= 1 + 0 (e~^) . 


(A25) 


(A26) 


Thus, for these low temperatures we recover as a good approximation the zero temperature expressions, and we obtain 
for the spinodal time a value very close to rj =0 . 


In summary, a good approximate expression for the spinodal time to order O 


l log log - 




IS 


T s = 


2 n 


log (t6) + i log 
i log (^f) + i log 


b» g (^) 


for (3 1 » /x , 
for (3~ l <€. g . 


(A27) 


This analytic expression is in good agreement with the numerical results, as can be seen in Table El 
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